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[X^ ■ Abstract 



We present a calculation of the next-to- leading order (0(q!^q;)) QCD corrections to heavy 
flavor photoproduction with longitudinally polarized beams. We apply our results to 
study the longitudinal spin asymmetry for the total charm quark production cross section 
which will be utilized by the forthcoming COMPASS experiment at CERN to obtain first 
direct information on the polarized gluon density A^f. We also briefly discuss the main 
theoretical uncertainties inherent in this calculation. In particular we demonstrate that 
the factorization scale dependence is considerably reduced in next-to-leading order. 



Despite significant experimental progress in spin-dependent deep inelastic scattering (DIS) 
the polarized parton densities A/ (/ = g, g, g) still remain considerably less known 
than their unpolarized counterparts /. In particular the polarized gluon density Ag 
is almost completely unconstrained in all recent next-to-leading order (NLO) analyses 
1^, ^, Q of presently available DIS data. This is due to the limited kinematical coverage 
of the data which does not allow for an unambiguous determination of Ag from scaling 
violations. An even more important role plays the lack of a momentum sum rule and the 
absence of any direct constraints from other processes like jet or direct photon production, 
both being available and important for the determination of unpolarized parton densities. 
Therefore a series of 'next generation' spin experiments will put a special emphasis on such 
exclusive measurements to provide further invaluable information for a more restrictive 
analysis of polarized parton distributions in the future. 

The first direct information on Ag is expected to be provided by the COMPASS fixed 
target experiment at CERN |^ from studies of the longitudinal spin asymmetry of total 
open charm photoproduction. Heavy quark {Q = c, b) production is generally considered 
to be one of the best options to pin down Ag since in leading order (LO) only the photon- 
gluon fusion (PGF) process 

19 ^QQ (1) 

contributes (an arrow denotes a longitudinally polarized particle). All studies of this 
process so far 0, 0) § were limited to LO estimates which, however, are known to be 
notoriously unreliable due to, e.g., their strong factorization/renormalization scale depen- 
dence. In addition, the already available unpolarized NLO corrections |10| turn out to 
be sizeable in certain kinematical regions and the 'clean picture' of Eq. ([^) is obscured by 
genuine NLO subprocesses with light quarks in the initial state. The knowledge of the 
polarized NLO cross section is thus mandatory for a meaningful extraction of A^f. It is 
the main purpose of this paper to provide for the first time the results of such a complete 
NLO {0{a'^a)) QCD calculation. Thereby we hope to provide a more reliable theoretical 
basis for the measurement of A^f by the COMPASS collaboration. 

The remainder of the paper is organized as follows: next we briefly sketch the most 
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important details of our calculation, mainly concentrating on complications which arise 
due to the presence of polarized particles in the initial state. We then apply our re- 
sults to present for the first time estimates for the longitudinal spin asymmetry of total 
open charm photoproduction in NLO QCD in the kinematical region accessible by the 
COMPASS experiment ||^. Finally, we address the relevance of the main theoretical un- 
certainties inherent in this calculation. For more details on our calculation, analytical 
NLO results, and studies of px and rapidity differential heavy quark cross sections, we 
refer the interested reader to Ref. [|l^ . 



The NLO QCD corrections to the PGF mechanism in Eq. ([l|) consist of three parts: 

(i) the one-loop virtual corrections, 

(ii) the real corrections with an additional gluon in the final state 

ig^QQg , (2) 

(iii) a new production mechanism appearing for the first time in NLO 

7qiq)^QQq{q) ■ (3) 



As is well known one encounters various types of singularities when calculating the loop- 
and 2^3 phase space integrals. Ultraviolet (UV) singularities which show up only in 
the virtual corrections (i) are removed by on-shell mass and wavefunction and (modified) 
MS coupling constant renormalization. In the latter case one also removes the gluon self- 
energy contribution due to the heavy quark loop in addition to the subtractions usually 
performed in the MS scheme (see, e.g., |lT], 0). This explicit decoupling of the heavy 
quark at low energy leads to a fixed flavor scheme with nif light flavors (active in the 
running of and in the parton evolution) and the produced heavy flavor. Infrared 
(IR) divergencies and double pole terms appearing in the non-abelian parts when IR and 
mass/collinear (M) singularities coincide cancel in the sum of (i) and (ii). Finally one 
is left with the M divergencies in (ii) and (iii) which are removed by the factorization 
procedure in the MS scheme. 
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One has to choose a consistent method of regularizing these singularities so that they 
become manifest. For this purpose we choose to work in the well-established framework 
of n-dimensional regularization. When calculating the required helicity dependent matrix 
elements we have to project onto the helicity states of the incoming particles. This is 
achieved by using the standard relations (see, e.g. 



(4) 



h ■ k2 

for incoming photons with momentum ki and helicity Ai (accordingly for gluons with k2 
and A2) and 

u{k2, h)u{k2, h) = ^ 1^2(1 - h-f^) (5) 

for incoming quarks with momentum k2 and helicity h (analogously for antiquarks). Using 
(^) and (|5|) we can calculate the contributions to heavy flavor photoproduction with 
unpolarized and polarized beams simultaneously by taking the sum or the difference of 
the helicity dependent squared matrix elements 

unpolarized : \M\'^ = \ (++) + (+-)] (6) 

polarized: A |M|' = ^ [|M|' (++) - |M|' (+-)] (7) 

where |M|^ [hi, h2) denotes the squared matrix element for any of the contributing sub- 
processes (|I]) and (i)-(iii) for definite helicities hi and /i2 of the two incoming particles. 
The possibility to obtain the unpolarized resultsQ 'for free' provides an important check 
on the correctness of our results. We fully agree with the unpolarized results presented 
in and [0. The abelian part of our result agrees analytically with the results of two 



recent calculations of the NLO QCD corrections to 77 QQ [|T4| as well. 

The presence of 75 and the totally anti-symmetric tensor e^yp„ in the polarized cal- 
culation due to @j and (^ introduces some complications because their purely four- 
dimensional origin allows no straightforward continuation to n 7^ 4 dimensions. We choose 
to work in the HVBM scheme [ITH] which was shown to provide an internally consistent 



"'^ Since each boson has n — 2 degrees of freedom, one has to average the spin of each incoming photon 
(gluon) with a factor l/(n — 2) replacing the 1/2 in (Q) for the unpolarized case. Furthermore, the 
inclusion of external ghost contributions (see, e.g., ||ll|) allows us to drop all terms other than g^i, in the 
symmetric part of (0). 
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extension of 75 and e^vpa to arbitrary dimension^. The price to pay is that apart from 
the usual n- dimensional scalar products k-p their respective (ra — 4)-dimensional subspace 
counterparts, usually denoted hy k ■ p ('hat momenta'), can also show up in |M|^ {hi, 
in the HVBM scheme. For the single-inclusive heavy quark production considered here. 



one can choose a convenient frame (see ITT] for details) where all non- vanishing (n — 4) 
scalar products can be expressed by a single hat momenta combination p^ = —p^. The 
latter quantity appears exclusively in the polarized 2 — 3 contributions and originates 
from scalar products of the momenta of the two not observed (integrated) final state 
particles in (^. However, terms proportional to p"^ still deserve special attention when 
performing the 2 — > 3 phase space integrations since the (n — 4)-dimensional subspace 
cannot be trivially integrated out as in an unpolarized calculation. 

The appropriately modified phase space formula for the present polarized calculation 
can be schematically expressed as follows 

dPSs = dPS3,.np(^l,^2) X / dx-= (8) 

with 

X = p / Pmax ~ 2~n • 2 n (^) 



where B is the Euler Beta function, m denotes the mass of the heavy quark, and S4 = 

s + ti + Ml is the sum of the three 2^2 Mandelstam variables. The angles 61^2 are 

introduced to parametrize the momenta of the two integrated partons (see [|TT[]). Since 

1 a;(«-6)/2 ( I ( I 

= { ^2 =^ A\r.'i.-2a c.\^2 n I ^2\i 5 (10) 



B{\,^)Jq I r \ (n-4)s2sin^^isin^^2/[4(s4 + m2)] 

Eq. (H) reduces to the well-known 'unpolarized' phase space dPSa^unp (see, e.g., |T^) for the 



vast majority of terms without a p^ dependence and for the rest, which is proportional to 
fp', introduces a.n ^ n + 2 shift in the angular integrals ([TT|) discussed below, guaranteeing 
collinear safety. Due to the extra powers of S4 all these additional integrals are infrared 
safe as well, implying that all contributions due to p^ are at least of (9(n — 4) (cf. Eq. (p!0|)) 
and hence drop out when the limit n 4 is takenQ. 

^This prescription is incorporated in TRACER which we have used for all trace calculations. 
•^This differs from a calculation involving only massless particles, where IR poles can occur. 
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The remaining phase space integration then proceeds as in the unpolarized case 
One has to make extensive use of relations between Mandelstam variables to reduce com- 
plex combinations to simpler ones by partial fractioning. In the end only integrals con- 
taining at most two angular dependent Mandelstam variables remain 

Jo Jo {a + bcos6i) (A -|- i? cos 6^1 + C sin cos 6*2) 
where a, b, A, B, and C are only functions of s, ti, ui, and m^. A useful list of integrals 
of the type ( pH]) can be found in [|17|, but we have checked all required formulae. The 



loop integrations have been performed using a Passarino-Veltman reduction |]T9| to scalar 



integrals which can also be found in |T^ and we have recalculated them as well. 

In the soft gluon limit for @ the 2^3 kinematics reduces to the usual 2^2 Born 
kinematics and soft gluon poles associated with I/S4 terms occur, since S4 — > when the 
final state gluon momentum becomes soft {k^ — > 0). To deal with these poles one can 
simply slice the phase space into a soft gluon (54 < A) and a hard gluon (54 > A) part 
0], where A has to be much smaller than s,ti,ui and m^. The soft gluon cross section, 



related to ds^, can then be obtained analytically and the singularities cancel in 



the sum with the virtual cross sections whereas the hard gluon part contains only the left 
over M singularities which are removed by factorization. 

As our final technical remark let us recall that the factorization 'counter cross section' 
for the light quark initiated subprocess (H) contains a part which can be schematically 
expressed as 

a.. 



+ ... 12 



where ® denotes a convolution, APg^ is the usual LO spin- dependent photon-to-quark 
splitting function and AFj^ represents the freedom in choosing a factorization prescription. 
In the MS scheme (which we use) AFg^ has the form 

AF^{x,f^f, /i) = APg,{x) (^^E - In 47r + In (13) 

where je is the Euler constant, fi is the mass scale introduced by n dimensional reg- 
ularization, and /i/ is the arbitrary factorization scale (in other schemes one subtracts 



different terms in (p!3|)). By applying (|T2|) one introduces the parton content of tlie real 
(on-shell) polarized photon |^ which is experimentally completely unknown for the time 
being. Strictly speaking a physically consistent (i.e., scheme independent) result for heavy 
quark photoproduction in can only be obtained for the sum of the 'direct' and 

'resolved' contribution due to the freedom in (|TB|). In the first case the photon acts as an 
elementary particle (as in the present calculation) whereas in the latter case it resolves 
into its parton content before the hard scattering takes place. The NLO corrections for 
the 'resolved' contribution are unknown so far for the polarized case and thus, for the 
time being, have to be estimated in LO (see below). 

Let us now turn to some numerical results and phenomenological aspects. We are here 
only interested in the experimentally most relevant total photoproduction cross section 
(for other results see [|ll|])- Let us start with the total photon-parton cross section in NLO 
which can be expressed in terms of scaling functions {i = g, q, q) 



A(Ji^(s, m^,/i/) 



^fl^iv) + 4vra. ( A/^«(r^) + Afi^{r,) In ^ 



(14) 



where A//^' and Afl^ , ^fl^ stand for the LO and NLO corrections, respectively, yU/ 
denotes the factorization scale, and t] = s/Aim? — 1. Note that for simplicity we have 
made the conventional choice /^r = A*/ in (HH) and in what follows (see for an 



independent variation of yUj and /i,.). The scaling functions can be further decomposed 
depending on the electric charge of the heavy and light quarks, cq and e^, respectively^: 

MaM = elAcg,{r]) (15) 
A/,7(^) = elAc,^{r]) + elAd^,{r]) (16) 

with corresponding expressions for the A/j^. 

In Fig. 1 we present Acg^ , Acf^ , and Ac^iy* as a function of the scaling parameter rj in 
the MS scheme (the threshold s = Aim? is located at 77 = 0). Also shown (dotted line) is 
the contribution to Acg^ from the hard gluon part (54 > A) aloneQ. From a comparison 

''Note that the interference term of the two possible production mechanisms for the subprocess (||) 
proportional to eqCq does not contribute to A/^-y in (p^, since it vanishes when integrated over the entire 
phase space as a consequence of the Furry theorem. 

^To define this quantity we follow the procedure in [p[ and explicitly add the A dependent terms 
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of this dotted curve with the dashed one showing the full NLO coefficient Ac^^ (i.e., 
the sum of the hard and soft plus virtual contributions), one can infer that the soft plus 
virtual contribution is almost negligible for f] ^ 1- For the interpretation of the results 
below it is important to notice that the LO coefficient (solid line) in Fig. 1 changes sign 
at ~ 3. Upon adding the NLO contributions, multiplied by a factor inag (see Eq. ([T^)), 
the zero is shifted towards ~ 1. We also note that for rj < 0.1 the 0{as) corrections 
dominate over the Born approximation when we include that factor. This enhancement 
near threshold is due to large logarithms from initial state gluon bremsstrahlung 
Furthermore in the threshold limit ?7 — the soft plus virtual piece tends to a constant 
('Coulomb singularity') while in LO Acf^^ 0. 

Fig. 2 compares in a similar way the light quark induced coefficients Acq^ , Acfcn j Ad^^ , 
and Ad^y as a function of t] in the MS scheme. Numerically they turn out to be much 
smaller than their gluonic counterparts shown in Fig. 1. Comparing our polarized and 
unpolarized results (for the latter see also Figs. 5, 7, and 8 in 0) for either the gluon or 
the quark coefficients, we find the same pattern: for — the unpolarized and polarized 
results become equal, e.g., Acg^ Cg^. A glance at Eqs. (||) and (|^) immediately implies 
that |Mj^|^ (H — ) 0. On the contrary, for asymptotically large energies r] oo the 
unpolarized NLO coefficients approach a large plateau value, except for rfq^ and d!^ , 
dominating over the LO result due to Feynman diagrams with a gluon exchange in the 
t-channel while all polarized coefficients tend to zero in that kinematical region. Thus 
one can infer that here |Mj^|^ (++) ~^ 1^17!^ (H — )• 

With the total partonic cross section in (0) at hand is is now straightforward to 

calculate the total hadronic heavy flavor photoproduction cross section via 

1 

Acr^p{Syp,m^,^f) = ^ J dxAarfy{xSyp,m^,fif)AfP{x,fij) (17) 

depending on the available photon-proton c .m.s. energy S^p. Of course an expression 
similar to (|l^) also holds for the unpolarized cross section cr^ with all polarized quanti- 
ties replaced by the corresponding unpolarized ones. Instead of measuring Ao"^ in ([T7|) 



(~ In A/m'^) of the soft plus virtual cross section to the hard part in order to cancel the dependence on 
the auxiliary quantity A numerically. 
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directly (which requires the determination of the absolute normalization), experiments 
will study the related longitudinal spin asymmetry defined by 



In Fig. 3 we show the LO and NLO predictions for the total charm asymmetry A^p, 
using m = 1.5 GeV, /i/ = 2m, and three different sets of polarized parton distributions 
P, ^ which mainly differ in Ag, in the ^/S^ region accessible by the upcoming COMPASS 
experiment at CERN Depending on the used muon beam energy (100 or 200 GeV) 
they will determine A^^ for one average value of \J S^p of about 10 GeV. The NLO 
corrections in Fig. 3 turn out to be sizeable, depend strongly on a/ S^p and do not cancel 
in the ratio (|l^) as one may naively expect. However, the physical origin of these large 
corrections is readily explained in terms of the LO and NLO coefficient functions (see 
Fig. 1 and 2). 

For \J S^p ~ 10 GeV one probes rj values from threshold up to 10. Note that in 
Acr^p, see Eq. (|17|), the coefficient functions for small values of rj (i.e., small s = xS^p) are 
convoluted with the gluon distribution at the smallest possible x {x > Aw? / S^p) where 
A(7 is usually large and vice versa. Since the NLO corrections to the gluonic coefficient 
functions, which dominate for small rj, decrease faster with increasing rj in the polarized 
case due to the zero, one should expect the NLO asymmetry to be somewhat smaller 
than the LO one in that particular S^p region if the LO and NLO /S.g are not too 
different. Moreover, the NLO shift of the zero towards smaller rj adds contributions with 
opposite sign in the convolution already for smaller S^p. This should also lower the NLO 
asymmetry with respect to the LO one for not too large S^p. Instead we find (see Fig. 3) 
that the NLO asymmetries for GRSV and GS (A) are larger than the LO ones for 
small S^p. But this is entirely due to the badly constrained A^f at large x and thus should 
not be taken too seriously: For S^p ~ 10 GeV the convolution (|T^ samples x > 0.1 and 
in both sets of polarized parton distributions p|, ^ the NLO Ag is much larger than the 
LO one (especially in the case of GS (A)). This is not the case for the unpolarized GRV 
distributions which we use to calculate cr^p in ([11). If we instead use NLO gluons 
for both the NLO and LO asymmetries, we find the expected behaviour discussed above 
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for GRSV and GS (A). This is demonstrated for the GRSV set by the dotted curve in 
Fig. 3. The situation for GS (C) 0, however, is even more complex since in this case 
the gluon oscillates as well in the covered x region. Also obviously the relative NLO 
corrections ('-ft'-factors') are always very large around zeros in the LO asymmetry, as for 
the GS (C) curve in Fig. 3 at around S^p ^ 12 GeV (similarly for GS (A) and GRSV 
at larger S^p ~ 30 GeV not shown in the figure), but this is natural for quantities that 
can change sign. 

Let us finally briefly discuss the importance of the main theoretical uncertainties in the 
calculation of Aa^^ or A^^ which may further complicate the extraction of /S.g. First of all 
Fig. 4 shows the dependence of Aa^p on the choice of the factorization scale in the relevant 
region /i^ = [m^, . . . , Aw?] for the GRSV standard set of polarized parton densities p[ for 
two different values of ■>/ S^p. The NLO results hardly depend on the precise choice of yU/ 
whereas the LO cross sections vary by as much as 50% in the shown range. This clearly 
underlines the importance (and usefulness) of the full NLO calculation. NLO results for 
Acr^p or Al^^ are thus much more trustworthy than previous LO estimates. 

As already mentioned above another complication in the extraction of /S.g arises due to 
the presence of the light quark induced subprocesses (§) which may serve as an important 
'background'. We have checked that for S^p ~ 10 GeV this contribution is fairly small 
(about 5%, except for the GS (C) set where the gluon contribution is close to zero at 
this \J S^p) and can possibly be neglected for a first determination of /S.g. However, 
with a better knowledge of the polarized quark densities in the future the light quark 
contribution should be subtracted. A similar remark holds for the importance of the 
possible 'background' from 'resolved' photons. Since the parton content of polarized 
photons is completely unknown so far, one has to impose some realistic models to 
estimate this contribution. In it was shown that even for rather large photonic densities 
this 'background' should also be very small in the 'COMPASS region'. Moreover, the 
unknown precise value of the charm quark mass m leads to shifts in Aa-yp and A'i^^ of 
about 30% when m is varied by 0.2 GeV around the central values of 1.5 GeV used in 
Figs. 3 and 4. 
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Finally, it should be stressed that for a meaningful and reliable extraction of Ag 
from a measurement of A'^^ (i.e., roughly a measurement of Ag/g) our knowledge of the 
unpolarized gluon density has to be improved as well. At \J S^p ^ 10 GeV one probes the 
gluon at X values larger than 0.1 where the uncertainty in (yf(x, /i^) is sizeable and non- 
negligible (see, e.g., P2| , In that particular x region g{x^ /i^) is mainly constrained by 
direct photon data, but with rather large theoretical uncertainties [^. Ideally, COMPASS 



should try to measure also the unpolarized charm photoproduction cross section, thereby 
reducing also our present ignorance on g at large x. Anyway, improvements seem to be 
mandatory for a reliable extraction of Ag. 

To summarize, we have presented the main results of a first complete NLO QCD 
calculation for heavy flavor photoproduction with longitudinally polarized beams. The 
NLO corrections are sizeable for the total longitudinal spin asymmetry in the energy 
range accessible by the forthcoming COMPASS experiment. But these corrections can 
be partly attributed to the rather different large x behaviour in LO and NLO of the 
presently available sets of polarized parton densities, which only reflects our complete 
ignorance of Ag. The theoretical uncertainties associated with the precise choice of a 
factorization scale were shown to be strongly reduced in NLO. The prospects for a first 
direct measurement of the polarized gluon density still seem to be promising but a careful 
and reliable extraction of Ag requires also an improved knowledge of the unpolarized 
gluon density at large x and of the precise value of the charm quark mass. Uncertainties 
due to the 'background' from 'resolved' photons or light quark induced processes should 
to be of less importance. 
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Figure Captions 



Fig. 1 The LO and NLO gluonic scaling functions in the MS scheme as a function of 
f] = s/Am? — 1 as defined in Eqs. (|1^) - (p^G]). 

Fig. 2 Same as in Fig. 1, but for the NLO hght quark coefficient functions. 

Fig. 3 The longitudinal spin asymmetry for total charm quark photopro duct ion in LO 



and NLO as defined in (|T8]) for m = L5GeV, jif = 2m and three different sets 
of LO and NLO polarized parton densities P, ^. The unpolarized cross section in 
( pISD was calculated using the GRV [^] densities. Also shown (dotted line) is the 
LO asymmetry using the NLO GRSV gluon distribution (see text). The vertical 
bar shows the estimated statistical uncertainty 5A'i^^ for such a measurement at 
COMPASS §. 

Fig. 4 The factorization scale dependence for the LO and NLO polarized total charm 
quark photoproduction cross section Acr^^ as defined in (|l^ using the LO and NLO 
GRSV standard densities 0, respectively, for m = 1.5 GeV. 
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Fig. 1 
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